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The discovery that the 11 pm and 12 pm window channels of AVHRR 
could be used to detect and even characterize the properties of cirrus (Inoue, 
1985 and others) stimulated the present study which re-examines the general 
multispectral approach for retrieving cirrus cloud top temperature and 
emissivity. The approach is based primarily on work done for METEOSAT 
by Szejwach (1982) and draws on a wealth of previous investigations of 
multispectral approaches (Coakley, 1983; Platt, 1983; Arking and Childs, 
1985). The generalized multispectral approach described here compliments 
the "C02 slicing method" used by Wylie (1988) and the bispectral threshold 
method used by Minnis et al (1988). 

The approach is based on the recognition that for a cloud system that 
is semi transparent in the infrared, the emission for channel i which senses 
radiation with wavelength %i is given by 

Ii = Isi + EcjAcdci ■ Isi) O) 

where Isi is the emission from the cloud-free region of the field of view and Ac 
is the fractional cloud cover within the field of view. Eci in (1) is an "effective 
mean" emissivity for the clouds within the field of view. Effects due to 
reflection of thermal radiation by the clouds are ignored in (1). Ici is the 
emission that the clouds would have were they blackbodies. Determining Ici 
is equivalent to determining the cloud top temperature. 

Consider now the behavior of observations at independent infrared 
wavelengths "hi and Xj under three conditions: 1) single-layered, opaque 

clouds, 2) single-layered, semitransparent clouds (circus) and 3) multi-level 
opaque clouds. In general, a cloud system will exhibit any combination of 
these three cases. Here, wavelengths that are independent are taken to be 
those that exhibit different sensitivities of the emission to temperature, i.e. 
are at substantially different wavelengths, or they exhibit different 
sensitivities to the dependencies of the optical properties to hydrometeor 
concentrations or both. 

For single-layered, opaque clouds, i.e. clouds that possess a single 
cloud top temperature, Is and Ici are constants at all wavelengths for the 
region containing the system. Likewise, because the clouds are opaque the 
emissivities, Ec, are also constants and equal to the maxi m u m values that the 
clouds can obtain. As a result, Ii is linearly proportional to I /. This linear 
relationship is shown in Figure 1 for the 11 pm and 12 pm radiances obtained 
from NOAA-9 on the afternoon of October 28, 1986 during the case study 
period. The data is for a 250 km region containing the array of lidars that 
participated in the case study (Sassen et al, 1989). The linear relationship 
that is shown in the figure indicates that the upper-level clouds were opaque 
at 11 and 12 pm at the time of the satellite overpass. On the basis of the 11 
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Figure 1. 11 (im and 12 (im AVHRR radiances for the afternoon NOAA-9 
pass on October 28, 1986. The data is for a 250 km region containing the 
array of lidars that participated in the Cirrus IFO case study. Each point in 
the figure represents a 4 km portion of the 250 km region. The linear 
relationship between the radiances indicates that from the view of the 
satellite, the clouds are opaque at 11 and 12 pm. The variability in the 11 
and 12 pm radiances indicate that the clouds are broken. 
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Figure 2. Same as Fig. 1 but for the 250 km region to the southwest of the 
lidar array. The curvature of the lower envelope of points indicates that the 
upper level clouds are semitransparent at some locations within the region. 
The intersections of the linear and nonlinear envelopes gives the cloud-free 
radiances, Is and the radiance associated with opaque clouds, Ic- The cloud 
top temperature could be estimated from Ic. 
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and 12 urn radiances, these clouds would not be interpreted as cirrus with 
missivities less than unity. While effects due to the reflection of theraml 
radiation have been ignored in this discussion, it is recognized that the 
maximum values of the emissivities attained by the clouds may be less than 
unity. The differential between unity and the maximum emissmty will 
affect the estimate of the cloud top temperature. 

For a single-layered, semitransparent cloud system Is and Id are 
constant for the region containing the clouds. The emissivities, however, can 
vary according to their dependence on hydrometeor concentrations. In 
general £ci is a nonlinear function of the hydrometeor concentration, and as a 
result £ci is a nonlinear function of £cj- Consequently, Ij will also be a 
nonlinear function of Lf . This nonlinear relationship is exhibited in Figure 2 
which shows 11 pm and 12 pm radiances from the October 28 NOAA-9 pass 
for the 250 km region southwest of the lidar array. The curvature exhibited 
by the envelope of points is due to the semitransparency of the upper level 
clouds. The linear relationship exhibited by the other envelope indicates that 
within this 250 km region, the upper-level clouds also become opaque at some 
locations. 

For clouds that are opaque but distributed in altitude, Ic will no longer 
be constant for the region containing the clouds. The emissivities, £c, will 
however be constant and near unity. As a result, the relationship between Ij 
and I/ will closely follow the dependence of the Planck function on 
temperature. That is, the relationship between If and Ij will again be 
nonlinear. In this case, however, the curvature will be opposite to that 
shown in Figure 2. 

One notes that taking the three cases together, one obtains a linear 
relationship between Ij and If for opaque, single-layered clouds, a nonlinear 
relationship between If and Ij for semitransparent, single-layered clouds and 
another nonlinear relationship for opaque, multi-layered clouds which follows 
the relationship given by the Planck function. The latter relationship may be 
calculated apriori. All three curves, however, intersect at Is ^hd *c where lc 
is that of the upper-level system which is semitransparent. So, in principle, 
one should be able to extract estimates for the cloud top temperature from 
the intersection. Having the cloud top temperature, one then estunates the 
distribution of e c A c in (1) from the distribution of observed intensities. 

While the results shown here were for 11 and 12 pm radiances, better 
definition of the cloud top temperature is probably obtainable using 3.7 pm 
radiances in combination with the 11 and 12 pm radiances. During the day, 
however, at least for the cases shown in Figs. 1 and 2, reflection of solar 
radiation at 3.7 pm by low-level water clouds makes the analysis untenable. 
At night, at least with the NOAA-9 AVHRR, instrument noise in the 3.7 pm 
channel also makes the analysis untenable. The identification of 
semitransparent systems using 3.7 pm radiances has, however, been noted 
elsewhere (Molnar and Coakley, 1985). 
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